size_genome=12071326
for length in {25,50,75}
do
for cover in {1,5,10,50}
do
for mean in {180,200,220}
do
cd "l=$length|f=$cover|m=$mean"
len_sam=`cat paired.sam |grep -v '^@'|wc -l`
m=`echo "scale=9; $len_sam*$length/$size_genome" | bc`
l=`echo "scale=9; e(-1*$m)" | bc -l`
c=`echo "scale=9; 1 - $l"|bc`
echo "$length,$cover,$mean,$len_sam,$m,$l,$c"
cd ../
done
done
done

